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I. INTRODUCTION 



The solutions of the non-relativistic and relativistic wave equations have attracted at- 
tentions in recent years. Furthermore, the analytical solutions of non-linear equations have 
also much attention. For example, the non-linear Landau-Lifshitz equation is solved to un- 
derstand the dynamics of Bose-Einstein condansates (BEC) in an optical lattice [1]. The 
two-component BEC with attractive interactions between atoms can be described by a non- 
linear Schrodinger equation (SE) called Gross-Pitaevskii equation [2]. 

Various methods are used to solve the Schrodinger equation based on perturbative and 
non-perturbative approaches such as the hypervirial-Pade summation method [3, 4], group- 
theoretical approach [5], Hill determinant method [6], and supersymmetric approaches [7]. 
The Klein-Gordon (KG) and Dirac equations are also studied for the Aharonov-Bohm (AB) 
potential [8], the AB plus the Dirac monopole potential [9], kink-like [10], Coulomb [11], 
vector plus scalar inversely linear potentials [12], PT-symmetric generalized Wood-Saxon 
(WS) [13], generalized Hulthen [14], and Rosen-Morse-type potentials [15]. These solutions 
are taken in general for constant mass [16,17]. On the other hand, position dependent mass 
case has also many application in different areas, such as impurities in crystals [18-20], the 
dependence of nuclear forces on the relative velocity of the two nucleons [21, 22], or the 
study of electronics properties of quantum wells and quantum dots [23], 3 He clusters [24], 
quantum liquids [25] and semiconductor heterostructures [26] . This is also to get the energy 
eigenvalues and eigenfunctions [27-31]. 

Here we intend to solve the KG-equation for the Woods-Saxon potential in the case of a 
exponentially mass distribution varying with position. In nuclear physics, the WS-potential 
is used to construct a shell model to describe the single-particle motion in a fusing system 
[32]. 

In the present work, we have obtained the energy spectrum and corresponding energy 
eigenfunctions by using the Nikiforov-Uvarov(NU)-method. We have also obtained the re- 
sults for the constant mass case. The NU-method is developed to solve the second order 
linear differential equations with special orthogonal functions. The method is based on 
solving the equation by reducing to a generalized equation hypergeometric type [40]. 

The organization of this work is as follows. In Section II, we solve the KG-equation in the 
case of the WS-potential for the spatially dependent mass by applying the NU-method, and 
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give the energy eigenvalues and the corresponding eigenfunctions. Our concluding remarks 
are given in Section III. 



II. NIKIFOROV-UVAROV METHOD AND CALCULATIONS 



A. 



Deformed Woods-Saxon Potential 



The KG-equation in one dimension for a particle reads 



dx 2 



1 



r 2 4 
[m c 



{E-Vf} 0(x) = O, 



(1) 



where E is the energy of the particle, c is the velocity of the light. The Woods-Saxon 
potential is 



is widely used in the coupled-channels calculations in heavy-ion physics. This model explains 
the single-particle motion during a heavy-ion collisions [32-35]. In this form of the potential, 
Vq is the potential depth, q is real parameter which determines the shape of the potential, 
and (3 is a short notation, i.e. (5 = 1/a, where a is diffuseness of the nuclear surface. 

Various mass-distributions are used in the literature. These are exponential, quadratic 
[28], inversely-quadratic [36], trigonometric mass-distributions [37], and mass function of the 
form m(r) = r a , is especially used for three-dimensional problems [37, 38]. Here, we prefer 
to use the following position dependent mass 



This provides us an exact solution of the effective KG-equation. tuq and mi in this distri- 
bution are two arbitrary positive parameters. m will correspond to the constant mass of 
the particle. The mass function is finite at infinity. 
Substituting Eqs. (3) and (2) into Eq. (1) we get 




(2) 




(3) 
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d 2 (j)(x) 
dx 2 



+ 



E 2 m 2 c 4 2EV /h 2 c 2 - 2m m 1 c 4 /ft 2 c 2 V 2 /h 2 c 2 - m 2 c i /n 2 c 2 



0(x) = O(4) 



h 2 c 2 h 2 c 2 (l + qe-P*) (l + qe-P*) 2 

To solve this equation, we use the transformation z — (1 + ?e~^ x ) -1 . By defining the 
following parameters 

we obtain 

d 2 <f){z) 1-2Z d<f){z) 1 r 2 2 2 21 , / x n / r n 

"dP" + 1^2 "dF + (J3P)2 [-^ -^-a 3 ]^) = o. (6) 

where Q 2 = l/(3 2 h 2 c 2 . Now to apply the NU-method [40], we rewrite Eq. (6) in the following 
form 

where <r(z) and ^(z) are polynomials with second-degree, at most, and f(z) is a polynomial 
with first-degree. We define a transformation for the total wave function as 

W=^(4 (8) 
Thus Eq. (7) is reduced to a hypergeometric type equation 

( t(^)^ ,/ (^) + T(z)ip'(z) + A^(z) = 0. (9) 
We also define the new eigenvalue for the Eq. (7) as 



A = A„ = -nr'- ^_il ff ",(n = 0,l,2,...) (10) 

Where 
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t(z) = f(z) + 2n(z). 



(11) 



The derivative of t(z) must be negative. A(A„) is obtained from a particular solution of the 
polynomial ip n (z) with the degree of n. ip n (z) is the hypergeometric type function whose 
solutions are given by [40] 



^ - mw y ' (z)p[z)l (12) 



where the weight function p(z) satisfies the equation 



±[v{z)p{z)]=T(z)p(z). (13) 



On the other hand, the function £(z) satisfies the relation 



?(z)/t(z)=*(z)/a(z). (14) 
Comparing Eq. (6) with Eq. (7), we have 

f(z) — 1 — 2z , a(z) = z(l — z) , a(z) = —a\z 2 — a\z — a\ (15) 

It becomes 



. , a'(z)—?(z) ,<?'{z) — t(z) .„ . , . . 

= 2 V ( 2 } " a(z) + ka(z) ' (16) 

or, explicitly 

ir(z) = T\J (af - k)z 2 + (a| + k)z + a§ , (17) 

The constant k is determined by imposing a condition such that the discriminant under 

the square root should be zero. The roots of k are k\^ = — 2a| T ^.a^A, where A = 
a/o| + a| + af. Substituting these values into Eq.(16), we get for tt(z) 
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ir(z) = ± 



(A — a 3 )z + a 3 , for k\ = —a\ — 1a\ + 2a 3 A 
(A + a 3 )z + a 3 , for k 2 = —a\ — 1a\ — 2a 3 A 



(18) 



Now we calculate the polynomial t(z) from tt(z) such that its derivative with respect to z 
must be negative, thus we take the first choice 



t(z) = 1 - 2(1 + A - a 3 )z - 2a 3 . 



(19) 



The constant A = k + tt'(z) becomes 



A = -a\ + (2a 3 - l)(A - a 3 ) , 



(20) 



and Eq. (10) gives us 



A„ = 2n(l + A - a 3 ) + n 2 - n . 



(21) 



Substituting the values of the parameters given by Eq. (5), and setting A = A„, one can find 
the energy eigenvalues as 



E n = Vo(Q 2 m omi c 4 - k/2)7 



=F 77y // 4/tm mic 4 - (k 2 /Q) 2 + 4m^(l/7 - Q 2 mfc 4 ) 



(22) 



where 



k= -(2n+l) 



i(2n + l) TV / l+4a? 



7' l=\l H + <>i- 



Q 2 v 2 + r 



■ (23) 



We see that the energy levels for particles and antiparticles are symmetric about 
v (Q m m,ic^-K/2) ^ rjnj ig g roun( j state energy is different from zero. To have a real energy 

spectra we impose 



k /2Q) + Q m m 1 c < K,m mic +m / / y, 



(24) 
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We plot four figures to present variation of first three energy eigenvalues as a functions 
potential parameters Vo and (3. Results are agreement with the ones obtained in the lit- 
erature [39]. M is the ratio mi/m , and 'p' and 'a' in the brackets represent 'particle'and 
'antiparticle' in figures. 

Now let us find the eigenfunctions. We first compute the weight function from Eqs. (15) 
and (19) 

p(z) = z- 2a *(l-z) 2A , (25) 

and the wave function becomes 

where b n is a normalization constant. The polynomial solutions can be written in terms of 
the Jacobi polynomials [41] 

^ n (z) = b n P^~ 2a ^(l-2z) 1 2A>-\, -2o 3 >-l. (27) 
On the other hand, the other part of the wave function is obtained from the Eq.(14) as 

£(*) = z a > (1 - z) A . (28) 
Thus, the total eigenfunctions take 

cf> n (z) = b' n (l- z) A z^ P (">-*«) (1 - 2z) (29) 
where b' n is the new normalization constant. It is obtained from 

f \K(z)\ 2 dz = l. (30) 
Jo 

To evaluate the integral, we use the following representation of the Jacobi polynomials [42] 
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p^) {z) = 



T(n + a + 1) T(ri + 1) Tjn + a + s + r + 1) z- 1 

n!r(n + (T + ^+ 1) ^ r(r + l)r(n-r + l) r(r + <j+l) ^ 2 ^' 



(31) 



Hence, from Eq. (30), and with the help of Eq. (31), we get 



F 2 J X Fl A s \b'J 2 f Z^+s (1 _ z) 2A dz = h (32) 

where -F^ 4 , and are two arbitrary functions of the parameters A, and a 3 , and given by 



?2A 



T(n + 2A + 1) 



n!r(n + 2A - 2a 3 + 1) ^ T(r + l)r(n - r + 1) 



r(n + l) T(n + 2A-2a 3 + r + 1) 



T(r + 2A + 1) 



(33) 



and 



i^-i^(n-m;r-s). (34) 
The integral in Eq. (32) can be evaluated from the definition of the Beta function [43] 



B{n,u) 
which gives us 



(35) 



Jo T(fl + i> + r + s) , 



(36) 



where fx = 2a 3 + 1, and v — 2A + 1. 

To get the energy eigenvalues for the constant mass case, we set q — 1, and m x = 



mi=0 



— =F 
2 



4V 2 + [v 7 /? 2 - 4V 2 - /3(2ra + l)] 2 



1 

16 



1/2 



(37) 



This is the same in Eq. (46) in Ref. [14]. 
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Since the wave function changes only with the parameter A, we simply get 



<f)^=°{z) = V' n (1 - z) A ' z a * pi** - 2z) , 
where the new parameter A' = \Ja\ — 2g 2 E/V — g 2 . 



(38) 



B. Non-PT Symmetric and non-Hermitian deformed Woods-Saxon Potential 

In this case, we take the potential parameters as Vo — > iVo, and (5 — > f3. So, the potential 
takes the form [44] 



V(x) = - 



iV 



1 + qe 



-fix 



(39) 



From Eq. (3), we obtain 



d 2 (f)(x) 

dx 2 



+ 



E 2 m 2 c 4 2iEV /H 2 c 2 - 2m m lC 4 /h 2 c 2 V 2 /h 2 c 2 + m 2 c 4 /h 2 c 2 



h 2 c 2 h 2 c 2 



+ 



(1 + qe-P x ) 



[1 + qe-?*) 2 



(f)(x) = . 



(40) 



By using the same coordinate transformation and defining the following parameters 



a 2 = Q 2 (E 2 -m 2 c 4 ) , -A 2 = Q (2iEV - 2m m 1 c ) , -A\ = Q 2 (-mfc 4 - V 2 ) , (41) 



we get 



Following the same procedure, we find the energy spectra 



(42) 



En 



iVo(2mom lC 4 - k"Q 2 ) 



2C 



^ ^V^(2m m 1 - k"Q 2 ) 2 - ((k" 2 - 4K,"Q 2 m omi c 4 + km 2 yQ 2 {Q 2 m 2 - k" - A 2 )) , 

(43) 



2C 
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where 



C = Q V + Q 2 m\c A ) , «" = - (2n + 1) 
and the corresponding total wave functions as 



1 (2n+l) + s]l + AA\ 



1 

+ 4 



(44) 



J n (z) = C (1 - z a3 P^ B ' ~ 2a ^ (1 - 2*) , 



(45) 



where B = \J A\ + A\ + of. It is seen that the energy eigenvalues are consist of the real and 
imaginary parts. For the constant mass case, we have 



E rn 1=0 = _ «V& ± 1_ ^ K „2Q2 {g 2 _ + + ^ , 



2 2Q 2 /c 
and the total eigenfunctions are 



(46) 



f,mi=0/ 



j) = 67 (1 - z a3 if 5 '- ~ 2a3 > (1 - 2z) , 



(47) 



where B' = a/* 3 ! ~ 2iEp 2 /V + p 2 . The energy spectra have real and imaginary part in the 
constant mass case. The imaginary part does not depend the quantum number n. The 
normalization constant b'" is also obtained in the same way. That is 



/~i2B (~i2B \yn\2 
Wr x ^ms \ u n I 



Jo 



2a3+r+s 



(1 - zf B dz = l, 



(48) 



The integral can be evaluated by using Eq. (35) 



1 Z 2a 3+r+s (1 _ Z) 2B dz = nm' + r + snn') 
y ' T(m> + n' + r + s) ' 

where m! = 2a 3 + 1, and n' = 2B + 1. Two functions G 2 f , and G 2 ® s are given by 



(49) 



G lr ^ F nr(^ A ^ 2B) , 

GZ - Fl A s {2A - 25) . (50) 
where the functions F 2 ; 4 , and are defined in Eqs. (33) and (34). 
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III. CONCLUSION 



We have solved the one dimensional effective mass KG-equation for the Woods-Saxon po- 
tential. The energy spectra and the corresponding wave functions are obtained by applying 
the NU-method. We have found a real energy spectra for the WS-potential in the position 
dependent mass case. To check our results, we have also calculated the energy eigenvalues of 
the particle and antiparticles for the constant mass limit. We have also studied the non- PT 
symmetric and non-Hermitian case. We have seen that the energy spectra have real and 
imaginary parts in this case. We have also obtained the energy spectra and corresponding 
eigenfunctions for the constant mass limit for this case. 
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FIG. 1: The dependence of first three excited states on Vq in the case of M = 0.04, and [3 = 0.1. 
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FIG. 3: The dependence of first three excited states on [3 in the case of M = 0.01, and Vq = 0.1. 
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FIG. 4: The dependence of first three excited states on (5 in the case of M = 0.01, and Vq = 1. 
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